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ABSTRACT 

We investigate the observed spectrum of cosmic-ray electrons and positrons from astrophysical 
sources, especially pulsars, and the physical processes for making the spectrum spiky or smooth via 
continuous and multiple electron/positron injections. We find that (1) the average electron spectrum 
predicted from nearby pulsars are consistent with PAMELA, Fermi and H.E.S.S. data. However, the 
ATIC /PPB-BETS peak around 500GeV is hard to produce by the sum of multiple pulsar contributions 
and requires a single (or a few) energetic pulsar(s). (2) A continuous injection produces a broad peak 
and a high energy tail above the peak, which can constrain the source duration (< 10 5 yr with the 
current data). (3) The H.E.S.S. data in the TeV range suggest that young sources with age less 
than ~ 6 x 10 4 yr are less energetic than ~ 10 48 crg. (4) We also expect a large dispersion in the 
TeV spectrum due to the small number of sources, that may cause the high energy cutoff inferred by 
H.E.S.S. and potentially provide a smoking-gun for the astrophysical origin. These spectral diagnostics 
can be refined in the near future by the CALET experiments to discriminate different astrophysical 
and dark matter origins. 

Subject headings: acceleration of particles - cosmic rays - pulsars: general 



1. INTRODUCTION 

Recently, the cosmic-ray positron fraction (the ratio 
of positrons to electrons plus positrons) has been mea- 
sured by PAMELA satellite (Adriani et al. 2008). The 
observed positron fraction rises in the energy range of 
lOGeV < e e ± < lOOGeV, contrary to the prediction 
of secondary positrons, which are generated from cos- 
mic rays propagating in the interstellar medium (ISM). 
The ATIC balloon experiment has also revealed that 
there is an excess above 300GeV and a possible peak 
at e e ± ~ 600GeV (Chang et al. 2008), which is also 
reported by PPB-BETS (Torii et al. 2008b). These ob- 
servations strongly indicate nearby sources of pairs 
within d ~ lkpc since high energy electrons/positrons 
lose their energy during propagation. Possible candi- 
dates include a pulsar (Shen 1970; Chi et al. 1996; Zhang 
& Cheng 2001; Grimani 2007; Kobayashi et al. 2004; 
Buesching et al. 2008; Hooper et al. 2009; Yuksel et al. 
2008; Profumo 2008; Malyshev et al. 2009; Grasso et al. 
2009), a microquasar (Heinz & Sunyaev 2002), a gamma- 
ray burst (GRB; Ioka 2008), a supernova remnant (SNR; 
Shen & Berkley 1968; Cowsik & Lee 1979; Erlykin & 
Wolfendale 2002; Pohl & Esposito 1998; Kobayashi et 
al. 2004; Shaviv et al. 2009; Fujita et al. 2009; Hu 
et al. 2009; Blasi 2009; Blasi & Serpico 2009; Mertsch 
& Sarkar 2009; Biermann et al. 2009) and dark matter 
annihilations/decays (Asano et al. 2007; Arkani-Hamed 
et al. 2009; Bergstrom et al. 2008; Hamaguchi et al. 
2008; Cirelli & Strumia 2008; Cholis et al. 2008a, 2008b; 
Chen et al. 2008, 2009a, 2009b; Chen & Takahashi 2008; 
Hisano et al. 2005, 2008a, 2008b, 2009; Ishiwata et al. 
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2008a, 2008b, 2008c; Zhang et al. 2008; March-Russell & 
West 2008; Hooper et al. 2008; Pohl 2009). Instead we 
might be observing the propagation effects (Delahaye et 
al. 2008; Cowsik & Burch 2009; Stawarz et al. 2009) or 
the proton contamination (Fazely et al. 2009; Schubnell 
2009). 

In order to discriminate different models of sources, 
an important diagnostic should be the spectral shape, in 
particular whether the ATIC/PPB-BETS peak is spiky 
or smooth. In this regard, it is remarkable that an astro- 
physical source can make a peak with a sharp cutoff that 
is similar to the dark matter predictions, if the source 
is a transient object like a GRB (Ioka 2008). However, 
other astrophysical sources like pulsars, SNRs or micro- 
quasars are not transient and expected to have a finite 
spread in the cutoff, as suggested by Ioka (2008). More 
importantly, due to the collimated emission, there are 
many off-axis pulsars that have not been observed via 
electromagnetic radiation, and we expect integrated con- 
tributions from multiple sources to the spectral shape, 
considering the birth rate of pulsars in our Galaxy. 

In addition, recently the Fermi Large Area Telescope 
has measured the electron spectrum up to ~ ITeV that is 
roughly proportional to ~ e^ 3 without any spectral peak 
as reported by ATIC/PPB-BETS (Abdo et al. 2009). 
The H.E.S.S. collaboration also provides the electron 
spectrum (Aharonian et al. 2008b, 2009), which is con- 
sistent with the Fermi result up to ~ ITeV and shows the 
steep drop of the flux above that energy. The Fermi data, 
however, should has a large systematic error in the high 
energy range (> 300GeV) where a significant fraction of 
electrons are removed to avoid a large hadron contami- 
nation, and so the real flux is estimated not by the pure 
experimental data but by the Monte Carlo simulations 
(Moiseev et al. 2007). On the other hand the ATIC data 
contains the larger statistical errors than the Fermi data. 
Therefore we cannot judge which observations are more 
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reliable so far. 

In this paper we investigate the effects of continu- 
ous and multiple pair injections on the observed elec- 
tron/positron spectrum. We show that the flux above 
the peak energy does not drop off abruptly but remain 
finite if the pair injection continues for a finite time and 
suggest that we may measure the source duration from 
the peak width. We also show that an average spectrum 
from multiple sources is relatively flat as reported by 
Fermi and the ATIC/PPB-BETS peak requires a single 
(or a few) extraordinary energetic source(s). We discuss 
the range of physical parameters of the sources (total 
electron/positron energy, the source duration, etc.) that 
are consistent with the current observational data. 

2. INJECTION MODELS AND CALCULATIONS 

2.1. Continuous e^ Injection from a Single Source 

We assume that a point-like source starts injecting e^ 
pairs at the time t = with total energy E e + ~ E e - at a 
distance d (~ lkpc) from the Earth. The observed elec- 
tron/positron spectrum after the propagation is obtained 
by solving the diffusion equation, 

- f = K(e e )V 2 f + g^{B(e e )f] + Q(t, r, e e ), (1) 

where f(t,r,e e ) is the distribution function of parti- 
cles at time t and position r with energy e e . Here 
K{e e ) — Kq(1 + e e /3GeV) 5 is the diffusion coefficient, 
B(e e ) is the energy loss rate, and Q is the injection 
rate of electrons/positrons. Hereafter we adopt Kq = 
5.8 x 10 28 cm 2 s _1 , 5 = 1/3 that is consistent with the 
boron/carbon ratio according to the latest GALPROP 
code, and B(e e ) = -be\ with b = lO^^GeV" 1 s" 1 which 
includes the energy loss due to synchrotron emission 
and inverse Compton scattering (Baltz & Edsjo 1999; 
Moskalenko & Strong 1998). 

Here we assume the continuous injection with a power- 
law spectrum: Q(t, r, s e ) oc Q (t)e~ a S(r — r ). We can 
obtain the observed spectrum for an arbitrary type of 
injection using the Green's function of the equation (|TJ) , 
derived in Atoyan et al. (1995), with respect to r and t: 
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x exp 
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where e e ,o — £e/[l — b(t — to)e e ] is the energy of elec- 
trons/positrons at the time to which are cooled down 
to e e at the time t, r = \r — ro|, and G — when 
e e fl is larger than the maximum energy of the injec- 
tion spectrum, e e!max . As B(e efi )/B(e e ) = (e e , /e e ) 2 < 
(e e ,max/£e) 2 < °o, there is no divergence in Eq. ([2]). We 
can approximate the diffusion length as 
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when e e ^ 3GeV and the diffusion coefficient is almost 
power-law K(e e ) ~ K (e e /3GeV) s . Here e cut = [b(t — 

to)}- 1 - 



Once we assume the injection rate Qo(to), we can ob- 
tain the observed electron/positron spectrum by inte- 
grating ([2]) over to- 



f(t,r,e e ) 



G(t,r,e e ;r, r )dr. 



(4) 



Note that the initial time of the integration should be 
set as U = max[0, t — 6 _1 (e~ 1 — e^^J]. 

We consider two types of continuous injection. One is 
the pulsar-type decay: 



Qo(t) oc 



(1 + r/ro) 2 



(5) 



This is the similar function of time as the spin-down 
luminosity of a pulsar with a surface magnetic field 



B = 8.6 x 10 1X P 
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where Piomsec is the pulsar period normalized by 10msec 
and Tq, 4 = r /10 4 yr (Shapiro & Teukolsky 1983). The 
other is the exponential decay: 



Qo(t) cx exp - 



Tn 4 
To 



(7) 



which may be realized by a pulsar that initially confines 
e ± in its nebula and releases them afterward, by a SNR 
that accelerates protons and continues to inject them to 
the surrounding dense gas cloud until it is destroyed, or 
by a microquasar ceasing its activity. In both types of 
injection, the characteristic time scale of the duration To 
is defined to be the time when the rate becomes four 
times smaller than the initial one. 

In Fig. 1 we show the electron plus positron flux re- 
sulting from above two injection models in addition to 
the transient model (to = 0) and the background 3 (dot- 
ted line). The remarkable point is that an astrophysical 
source can make a spectral peak that is similar to the 
ATIC/PPB-BETS excess and also to the dark matter 
case (Ioka 2008). The peak energy is determined by the 
age of the source £ age as 



£e,pcak 



1 



(8) 



because the electrons/positrons with initially higher en- 
ergy cool down via synchrotron emission and inverse 
Compton scattering within time i agc . We can inversely 
estimate the source age as i ag0 ~ 5 x 10 5 years from the 
peak energy for £ e ,max ^ ITeV. Note that the peak flux 
is almost independent of the distance r if it is smaller 
than the diffusion length (~ lkpc in our case). 

As is clear from Fig. 1, the spectral cutoff becomes 
shallower for the continuous injection models than the 
transient one (r = 0; short dot-dashed line). This is 
because the significant fraction of e pairs are produced 
recently (i.e. injected long after the birth of the source) 
and they have shorter time for the energy loss via syn- 
chrotron emission and inverse Compton scattering. Then 

3 For the background shown in the following plots, we adopt the 
fitting functions in Baltz & Edsjo (1999) by reducing the primary 
e~ flux, which is conventionally attributed to supernova remnants, 
by 30% because the fitting functions provide larger flux than the 
ATIC data even without other contributions. 
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Fig. 1. — The electron plus positron flux predicted from a source 
that continuously injects pairs for a finite duration to = 10 5 years 
with the exponential decay in Eq. (7) (thin solid line), and its sum 
(thick solid line) with the background (dotted line), compared with 
the ATIC/PPB-BETS/H.E.S.S. /Fermi data. We also show the 
pulsar-type injection in Eq. (5) with to = 10 5 years (long dashed 
line) and to = 10 4 years (double dashed line), in addition to the 
transient injection (to = 0; short dot-dashed line). We assume that 
a source at r = lkpc from the Earth a time t age = 5.6 X 10 5 years 
ago produces e pairs with total energy E e + = E e - = 0.8 X 
10 5( ^erg and spectral index a = 1.7 up to e e .max — lOTeV. 

their energy is still higher than the peak energy when 
they reach the Earth, and they produce a broader peak. 

The thick solid line represents the total (the primary 
plus background electron and positron) flux assuming 
that the source starts emitting e pairs with total energy 
~ 10 50 erg, a power-law index a ~ 1.7 and a maximum 
energy ~ 5TeV at a distance ~ lkpc from the Earth a 
time i ag0 ~ 5x 10 5 yr ago, and decays exponentially with 
the duration of To ~ 10 5 year. This model looks better for 
the ATIC/PPB-BETS peak, though we cannot conclude 
that the duration is finite with the current data. The 
positron fraction predicted from this parameter set is also 
consistent with the PAMELA results, in almost the same 
way as Fig. 1 of Ioka (2008). 

In the case of pulsar-type injection, there is another 
interesting spectral feature resulting from a long dura- 
tion. In Fig. 1, the high energy tail above the peak 
energy is more enhanced for the long duration case 
(to = 10 5 years, double dashed line) than the short du- 
ration case (to = 10 4 years, long dashed line). This is 
because the longer the duration of injection is, the larger 
fraction of e ± pairs are freshly produced and they do 
not lose their energy during the propagation so much 
(see also Atoyan et al. 1995). Especially, the flux of the 
long duration model may exceed the H.E.S.S. observa- 
tions around ~ 4TeV if we add the background (dotted 
line) while that of the short duration model does not. 
As the errorbars are still large, however, we should await 
future observations. 

2.2. Multiple e Injections: Average Flux and Its 
Dispersion 

Next, let us consider multiple sources. We expect sev- 
eral younger or older pulsars than that in Fig. 1 of age 
t agc ~ 5 x 10 5 yr, considering the local birth rate of pulsars 
~ 10~ 5 yr _1 kpc~ 2 (Narayan 1984; Lorimer et al. 1993). 
Moreover, the total energy E e + + E e - in Fig. 1 is as 
large as the rotation energy of a pulsar E rot with a pe- 
riod of ~ 10msec. This is comparable with the fastest 



initial spin estimated from the observations of radio pul- 
sars (Kaspi & Helfand 2002), so the pair output efficiency 
f e = (E e + + E e -)/E Tot may be too large ~ 100% to ac- 
count for the excess with a single pulsar. 

We can calculate the average electron and positron 
spectrum by considering the nearby multiple pulsars with 
a certain birth rate in the following way. Once e e is fixed, 
we can neglect the contribution from pulsars older than 
~ l/(be e ) and farther than ~ e?diff ~ 2*JK(e e )t as is 
obvious from the functional form of Eq.((2]). Then the 
average flux can be calculated by: 



/ave(^e) : 



l/(b6e) 



dt 



QoR 



2irrdrf(t, r, e)R 



= N(e e ) x /i,avc(e e ), 
where R is the local pulsar birth rate (yr^ 1 kpc -2 ), 

~l/(be e ) 



(9) 
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dt 



2irK(e e )R 
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(be e y 



(10) 



is the number of pulsars which contribute to the flux at 
the energy e e , and /i, a vc(e e ) = /avc(e e )/A(e e ) is the av- 
erage electron flux per pulsar. Here we adopt the value of 
R as the birth rate per unit surface area because pulsars 
are born from a disk whose thickness (~ 200 — 300pc) 
is much smaller than the diffusion length of cosmic-rays 
(~2-3kpc). 

Fig. 2 and Fig. 3 show the average electron spectra and 
the positron fraction, respectively, obtained by assuming 
that each pulsar emits electrons with the total amount 
of energy of ~ 1 x 10 48 erg, the spectral index of ~ 1.9, 
and the birth rate of R ~ 1/(1.5 x 10 5 )yr _1 kpc -2 (thick 
solid lines). 

In addition, we can calculate the dispersion of the num- 
ber of pulsars AN(e e ) from the average N(e e ) for each 
energy bin as 



AN(e e ) ~ y/N(e e ), 



(11) 



which is based on the Poisson distribution of nearby pul- 
sars. Then we can estimate the flux dispersion as 



A/ avc (£ e ) ~ /l,avc(e e )v / ^(£e) = /ave(ee)/ ^N(e e %12) 

From Fig. 2 and 3 we can see that the average spec- 
tra are basically consistent with Fermi, H.E.S.S. and 
PAMELA data. In the high energy range (e e > TeV), 
the dispersion from the average flux become significant. 
This can be interpreted as follows. The pulsars which 
contribute to the electron and positron flux in such a 
high energy band should be young, 



iage „ -rj 



3.1 x 10 5 yre 



(— y 1 

\TcVJ 



and close to the Earth, 



< 2^K{e e )t agc ~ 1.3kpc (^) 



-1/3 



(13) 
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where we adopt S = 1/3. The number of such pulsars 
should be as small as 



N(e e ) 



VTeV/ V 1/(1- 



R 



/(1.5 x 10 5 )yr- 1 kpc" 



(15) 



Therefore, in the TeV range few pulsars can contribute 
to the electron/positron flux. This small number of pul- 
sars may naturally account for the spectral cutoff around 
~ TeV energy, which has been inferred by the H.E.S.S. 
observations. Strictly speaking, this estimation of the 
pulsar number dispersion is not correct in the energy 
range of e e > 3TeV, where N(e e ) < 1 and the statistical 
arguments become meaningless. However this interpre- 
tation of the spectral drop around this energy is still 
qualitatively correct. 

Moreover, Fig.2 shows that the ATIC/PPB-BETS 
peak flux [e e ~ 545GeV) is much larger than the av- 
erage flux added with the dispersion flux A/ avo at the 
same energy bin. In fact, the separation between the av- 
erage flux and the ATIC data of the peak flux at that 
energy is ~ 10A/ avc . Then, if all pulsars emit electrons 
with the total energy of ~ 10 48 erg, the number of pulsars 
which contribute to the energy bin of the ATIC/PPB- 
BETS should be unrealistically large at the 10a level. 
This means that if the ATIC/PPB-BETS peak is real, it 
does not seem to be produced by the collective contri- 
bution from multiple pulsars with the moderate amount 
of electron energy (~ 10 48 erg) but by a single (or a few) 
energetic pulsar(s) (~ 10 49_50 erg) 4 . 

The discussion above is about the fluctuation of the 
number of pulsars with a certain birth rate, and it is 
based on the Poisson statistics. Strictly speaking, in or- 
der to discuss the cosmic ray electron/positron fluctua- 
tions due to the random injections, one should evaluate 
not the dispersion of the source number but the disper- 
sion of the electron/positron flux at each energy bin as 



A/, 
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avc 



l/(tee 
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2vrdrf 2 R-N(e e )f 2 (l&) 



The first integral in Eq. (|16p contains, however, a se- 
rious divergence because of the large (but improbable) 
contribution from very young and nearby sources (Lee 
1979; Berezinskii et al. 1990; Lagutin and Nikulin 1995; 
Ptuskin et al. 2006). In order to obtain the realistic es- 
timate of the flux dispersion, we introduce a lower cutoff 
parameter r c to the time integral. Then we have 
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l6ir 2 K(e e ) 2 T c 



(17) 



Following Ptuskin et al. (2006), we adopt the cutoff 
parameter as 



: [AnRK(e e )] 



-1/2 



: 10 5 yr 



( £e ) 
V500GeV/ 



-1/6 



R 



,500GeV7 ^1/(1.5 x lO^yr-ikpc' 2 
which takes into account the absence of very young and 

4 Since the multiple contributions tend to make the spectrum 
softer, it is possible to fit the ATIC/PPB-BETS spectrum with 
multiple pulsars by using the harder spectral index a accordingly. 
However, in order to fit the PAMELA spectrum in the lower energy 
range at the same time, such a hard spectrum is not favored. 
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Fig. 2. — The average electron plus positron flux (thick solid 
line), the flux with the standard deviation (/ ave ± A/ aV e, up- 
per/lower thick dashed lines) with the background (dotted line) 
predicted from the local pulsar birth rate of ~ 1/(1.5 X 10 5 ) per 
year per kpc 2 and the total electron energy of ~ 10 48 erg (thin 
solid line and upper/lower thin dashed lines for the same set of the 
spectra without the background), compared with the ATIC/PPB- 
BETS/H.E.S.S. /Fermi data. The flux from a young source (~ 
6 X 10 4 years) with the energy of ~ 2 X 10 48 erg (long-dashed line) is 
also shown. We assume that each source emits electrons/positrons 
with a power-law index a = 1.9 up to e e ,max = lOTeV. 



nearby sources. This choice of r is reasonable as long 
as this time is much shorter than (6e e ) _1 , which means 
e e <C 4.4TeV. Then the ratio of the flux dispersion to 
the average flux ((9]) can be expressed as 



A/a 



^0.21 



500GeV 
R 



5/12 



1/(1.5 x 10 5 )yr- 1 kpc" 



-1/4 



(19) 



We can see that these two 'dispersions' (Eq. (fT2| and 
Eq.(|19p) give the similar results around the energy of the 
ATIC/PPB-BETS peak. Therefore, from either of the 
above discussions, we can say that the ATIC/PPB-BETS 
data of the spectral peak is so largely separated from the 
average flux that they do not seem to be produced by the 
multiple contribution from nearby pulsars with moderate 
energy. 

We should note that the H.E.S.S. data put constraints 
on the total e pair energy from young sources. We plot 
in Fig. 2 the electron spectrum from the source with the 
age of ~ 6 x 10 4 years so as not to exceed the observational 
upper limit inferred by the H.E.S.S. data in the TeV 
range (long-dashed line). We find that the total energy 
of such young sources should be, if exist, < 2 x 10 48 erg 
which is two orders-of-magnitude smaller than the energy 
of the source making the ATIC/PPB-BETS peak. 

3. DISCUSSION AND CONCLUSION 

We investigate the astrophysical origin for the 
PAMELA and ATIC/PPB-BETS excesses and in par- 
ticular the effects of the finite duration and the multiple 
sources on the electron and positron spectra, as expected 
for pulsars, SNRs and microquasars. We find the follow- 
ings: 

(1) A non-transient source can make a spectral peak 
that is similar to the ATIC/PPB-BETS excess (see 
Fig. 1) around the peak energy in Eq. (8J. The peak 
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Fig. 3. — The total positron fraction resulting from the average 
spectrum (thick solid line) and the dispersion (thick dotted lines), 
which have the same parameters as in Fig. 2, and the background 
(dotted line), compared with the PAMELA data. Note that the 
solar modulation is important below ~ lOGeV. 



is generally broad with a width 



Ae. 



e.pcak 



^e,peak 



~ 10% 7-0,4 ( 



i' age 



\ 10 5 years 



(20) 



which could provide a method to measure the source du- 
ration To by the Fermi satellite (an energy resolution of 
5-20% in 20GeV-lTeV range; Moiseev et al. 2007) or 
the future CALET experiments (a few % above lOOGeV; 
Torii et al. 2008a). Although Atoyan et al. (1995) have 
already pointed out the effects of finite duration of the 
source on the electron spectrum, they only mention the 
enhancement of the high energy tail above the spectral 
peak (see below) and never discuss the peak width. Note 
that the peak width is also produced by the spatial fluc- 
tuation of Galactic magnetic field and the photon den- 
sity because the energy loss rate of fluctuates during 
the propagation, as estimated in Ioka (2008) (see also 
Malyshev et al. 2009). We also note that the peak be- 
comes smoother if the injection rises gradually in its ini- 
tial stage. 

(2) The spectrum from a long duration source has a 
high energy tail above the peak energy (see Fig. 1). Espe- 
cially the flux of this tail plus the background may exceed 
the H.E.S.S. data points when assuming a pulsar-type de- 
cay with a duration To > 10 5 years. This implies that the 
source making the ATIC/PPB-BETS peak is not likely 
a single pulsar with magnetic fields weaker than a few 
times 10 11 G. The existence of this tail has been already 
pointed out before (Atoyan et al. 1995). However, we 
firstly present the quantitative argument for the obser- 
vational limit of the duration of the electron/positron 
source in the context of the high energy tail thanks to 
the observational developments in the TeV range. One 
should note that we cannot rule out the long-duration 
pulsar model if the maximum energy of injected e^- pairs 
is smaller than < TeV, or the injection is not the pulsar- 
type in Eq. (5) but the exponential-type in Eq. (7), for 
example. The latter is possible if high energy pairs gen- 
erated in the pulsar magnetosphere are not injected into 
the space instantaneously but initially confined in a pul- 
sar wind nebula (Chi et al. 1996) and they diffuse out 
after the nebula gets broken. 

(3) The H.E.S.S. data suggest that young sources with 
age less than 6 x 10 4 yr should be, if exist, two orders- 



of-magnitude less energetic than the source making the 
ATIC/PPB-BETS peak. Note that the lifetime of the 
pulsar nebula is around ~ 10 5 yr and younger pulsars may 
not be able to contribute by the cosmic-ray confinement 
in the nebula. 

(4) The average electron spectrum and positron frac- 
tion is well consistent with the H.E.S.S. /Fermi and 
PAMELA data, respectively, taking into account the 
dispersion predicted from the total electron energy per 
pulsar of ~ 10 48 crg with the local birth rate of ~ 
l/10 5 yrs/kpc 2 . Especially, when e e > TeV, we expect 
a large dispersion of the electron flux because of the 
small number of sources which are young and close to 
the Earth and can significantly contribute to that energy 
range. This fact can naturally account for the spectral 
drop around > TeV indicated by the H.E.S.S. observa- 
tions. Note that the value of the total electron energy per 
pulsar adopted here is within reasonable range. In fact, 
the pulsar whose initial spin period is around ~ 10msec 
can emit electrons and positrons with energy ~ 10 48 erg 
if we assume the efficiency of f e ~ 1%, which seems to 
be reasonable (Hooper et al. 2009). 

Moreover, we show that the ATIC/PPB-BETS data 
point showing the peak at e e ~ 545GcV is largely sep- 
arated from the average flux, when considering the the- 
oretical dispersion from the average. This fact suggests 
that the peak is hard to produce by multiple contribu- 
tions and requires a single (or a few) extraordinary pul- 
sar(s) whose total electron/positron energy is about a 
hundred (several tens) times larger than that of ordinary 
pulsars. Here we estimate the dispersion of the elec- 
tron/positron flux based on the analytical expressions 
(Eq. (9), (12) or (17)) using the averaged local birth rate 
of pulsars. This method enables us to take into account 
the off-axis pulsars whose existence is suggested by the 
observed pulse shape of pulsars, and it is different from 
the method used in Malyshev et al. (2009) who calcu- 
late some realizations of the spectra predicted from the 
known pulsars in the ATNF catalogue. 

Note that the different choice of the diffusion coeffi- 
cient K(s e ) would change the results quantitatively. The 
smaller K makes the diffusion length r<jiff smaller, and 
the particle density inside that radius gets higher, being 
proportional to rj ; ^ (see Eqs.(2) and (3)). For different 
K instead of K, we can apply our results by re-scaling 
the distance of each pulsar and the total injection 

energy as d -► d^K/K and E tot -> E tot (K /K ) 3 / 2 , re- 
spectively. 

In our calculations we evaluate the dispersion of the 
electron flux due to the random birth of nearby pulsars in 
time and space having uniform total energy and injection 
index. In the case that these pulsars have a distribution 
of energy with a dispersion of 5E, the total dispersion 
of the energy is averaged as ~ 5Ey/N(e), and when the 
electron energy is smaller than ~ TeV (i.e. N(e) is much 
larger than unity) the total dispersion is suppressed com- 
pared to the total flux N(s e )fi, ave . The spectral index 
of the injected electrons should also be varied. However, 
the dispersion of the flux is almost determined by the 
amount of the electron energy emitted from pulsars, and 
the fluctuation of the index would not contribute to the 
flux dispersion so much. 
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The spatial variation of the energy loss rate and the dif- 
fusion coefficient can also affect to the observed electron 
flux or positron fraction. The energy loss rate b can fluc- 
tuate along the propagation path of electrons because of 
the inhomogeneities of the radiation and magnetic field, 
and then the cutoff shape of the resulting electron spec- 
tra would be broadened according to the amplitude of the 
fluctuation. Such a feature may be resolved by the future 
CALET experiment (see Ioka 2008). On the other hand, 
the effects of the spatial variation of the diffusion coef- 
ficient are considered in Cowsik & Burch (2009) in the 
context of "Nested leaky box model". In this model the 
positron fraction can be explained as a result of the differ- 
ent diffusion coefficient between the source-surrounding 
region and the general interstellar space. 

We can expect gamma-ray emission from high energy 



pairs. Especially, the number of such energetic ob- 
jects can be simply estimated as ~ (10kpc/lkpc) 2 = 100. 
This is comparable with that of TcV unidentified sources, 
which have no clear counterpart at other wavelengths 
(Aharonian et al. 2005, 2008a; Mukherjee and Halpern 
2005; Ioka & Meszaros 2009), implying some connections 
between them. 
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Grant-in- Aid for Science Research, Japan Society for the 
Promotion of Science (No. 18740147 and No. 19047004 
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